## Examples of the Möbius mu function

from PyM import *

m = mu_moebius

show(m(1), m(101), m(35), m(101101101))

show(ifactor(101101101))

# Example of Möbius inversion formula
n = 9*5*7

D = divisors(n)

show(phi_euler(n))

S = sum([m(d)*n//d for d in D])

show(S)